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The IceCube collaboration has recently reported the observation of two events with energies in 
excess of 1 PeV. While an atmospheric origin of these events cannot be ruled out at this time, this 
pair of showers may potentially represent the first observation of high-energy astrophysical neutrinos. 
In this paper, we argue that if these events are neutrino-induced, then the neutrinos are very likely 
to have been produced via photo-meson interactions taking place in the same class of astrophysical 
objects that are responsible for the acceleration of the ~10 17 eV cosmic ray spectrum. Among 
the proposed sources of such cosmic rays, gamma-ray bursts stand out as particularly capable of 
generating PeV neutrinos at the level implied by IceCube's two events. In contrast, the radiation 
fields in typical active galactic nuclei models are likely dominated by lower energy (UV) photons, and 
thus feature higher energy thresholds for pion production, leading to neutrino spectra which peak at 
EeV rather than PeV energies (models with significant densities of x-ray emission, however, could 
evade this problem). Cosmogenic neutrinos generated from the propagation of ultra-high energy 
cosmic rays similarly peak at energies that are much higher than those of the events reported by 
IceCube. 
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I. INTRODUCTION 

Very recently, the IceCube collaboration announced 
their observation of two events which could potentially 
represent the first detection of high-energy astrophysical 
neutrinos [T]. The analysis under consideration was de- 
signed to search for very high-energy (> PeV) and high- 
quality shower events over a period of 672.7 live days 
between 2010 and 2012. Over this time, the IceCube 
experiment was in its nearly complete (79 string) and 
complete (86 string) configurations. The two showers in 
question were observed on August 9, 2011 and January 3, 
2012, and each have an energy of approximately 1-2 PeV. 
Both showers were fully contained within the volume of 
the detector, and there are no indications of any instru- 
mental problems, or of any connection with atmospheric 
muons. 

At energies below ~1 PeV, atmospheric neutrinos con- 
stitute IceCube's primary background. Conventional at- 
mospheric neutrinos, which are produced in the decays 
of pions and kaons present in cosmic ray induced cas- 
cades in the Earth's atmosphere, follow a spectrum which 
falls-off rapidly with energy (following a power-law spec- 
trum with an index of approximately 3.7). Although this 
background is negligible at energies above several PeV, 
it cannot be entirely discounted in the case of the anal- 
ysis at hand. In particular, the IceCube collaboration 
expects 0.14 background events (including those from 
conventional atmospheric neutrinos and from misiden- 
tified atmospheric muons) over the time period covered 
by their analysis. The observation of two events with 
an expected background of 0.14 constitutes a P- value of 
0.0094, corresponding to a significance of 2.36er. 

In addition to the so-called conventional atmospheric 



neutrinos, a somewhat harder spectrum of prompt atmo- 
spheric neutrinos produced in the decays of heavy flavor 
(charm, bottom) mesons is also expected. Estimates of 
the event rate from this component of the atmospheric 
neutrino spectrum are less certain, although the lack of 
any correlation between the two observed events and hits 
in the IceTop surface array appears to disfavor this in- 
terpretation (efforts to quantify IceTop's veto efficiency 
are currently ongoing). 

Given the quite modest statistical significance (2.36c) 
represented by these two events, we cannot at this time 
be confident that they are of cosmic origin. That being 
said, if they are in fact cosmic neutrinos, then a few more 
years of data should reveal several more PeV showers. At 
a continued rate of 2 events per 672.7 live days, we es- 
timate that this signal will become inconsistent with an 
atmospheric origin at the 5tr level after a total of approx- 
imately 7 years of data taking. Searches for muon-track 
events or partially contained showers could also shed a 
great deal of light on the question of whether these events 
in fact originate from beyond the Earth's atmosphere. 

In this paper, we assume that this pair of intriguing 
events are astrophysical in nature and discuss the types of 
astrophysical sources from which they could potentially 
originate. In Sec. [TTJ, we discuss in general terms the as- 
trophysical production of PeV neutrinos and show that 
these two events - if, in fact, astrophysical in nature - 
are likely to originate from the same class of objects that 
produce the bulk of the 10 17 eV-scale cosmic rays. In 
particular, for such sources to produce the PeV neutrino 
flux implied by IceCube's two events, it is only necessary 
that on the order of ten percent of the energy in ~10 17 eV 
cosmic rays is lost to photo-meson interactions. In this 
context, gamma-ray bursts are a particularly attractive 
possibility, which we discuss in Sec. |III| We consider other 
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possible sources such as active galactic nuclei in Sec. |IV| 
but find these scenarios less compelling unless the objects 
in question are surrounded by significant densities of high 
energy (x-ray) radiation. We also present arguments for 
why these neutrinos are unlikely to be produced in the 
propagation of cosmic rays (cosmogenic neutrinos) or re- 
sult from the decays of ultra-high energy neutrons. In 
Sec. |Vj we make predictions for future observations at 
IceCube and, in Sec. VI we summarize our results and 
conclusions. 



II. THE ASTROPHYSICAL PRODUCTION OF 
PEV NEUTRINOS 



PeV-scale astrophysical neutrinos can be produced 
through three primary processes. Firstly, proton colli- 
sions with energetic photons can generate charged pions 
(photo- meson production), which yield neutrinos in their 
decays. To exceed the threshold for pion production, 
however, one must consider circumstances in which the 
target radiation is fairly energetic (10's of eV or above). 
Secondly, collisions between energetic protons (or neu- 
trons) and gas can easily exceed the threshold for pion 
production. In sources of protons with enough energy 
to yield PeV neutrinos, however, the number density of 
sufficiently energetic photon targets present will almost 
certainly be much larger than that of nucleons. Thus 
we only expect pp collisions to dominate neutrino pro- 
duction in sources of lower energy neutrinos. And lastly, 
PeV electron anti-neutrinos can result from the decays of 
ultra high-energy neutrons. As only a small fraction of 
a neutron's energy goes into its neutrino decay product, 
however, neutron decay anti-neutrinos are not predicted 
to be produced in sufficient numbers to account for the 
two events recently reported by the IceCube collabora- 
tion (see, for example, Refs. [2JI3]). 

With these considerations in mind, we expect cosmic 
ray collisions with energetic photons to be the dominant 
mechanism behind IceCube's two reported events (as- 
suming they are in fact astrophysical and not terrestrial 
in nature). And while one could also contemplate more 
exotic sources, such as the decays of ultra-heavy parti- 
cles, we confine ourselves to conventional astrophysical 
scenarios in this paper. 

Each of the three production mechanisms described 
above involve very high-energy cosmic rays, providing 
us with a direct connection between the observation of 
high-energy neutrinos and the observed cosmic ray spec- 
trum [4] . More specifically, one can place an upper bound 
on the diffuse neutrino flux that results from cosmo- 
logically distributed, optically-thin cosmic ray acceler- 
ators. This argument, originally presented by Waxman 
and Bahcall in connection with cosmic rays of the highest 
energies [5] , begins with the (cosmologically) local energy 
injection rate of ultra high-energy (10 19 — 10 21 eV) cosmic 
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where an injected energy spectrum cx E 2 has been 
assumed. The energy density of neutrinos produced 
through the photo-meson interactions of these protons 
can be directly tied to the injection rate of cosmic rays: 
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where tn is the Hubble time and e T is the fraction of the 
energy which is injected in protons lost to photo-meson 
interactions. The factor of 3/8 comes from the fact that, 
near the threshold for pion production, roughly half the 
pions produced in photo-meson interactions are neutral 
and do not generate neutrinos, and three quarters of the 
energy of charged pion decays (ir + —> ^jl + v^ — > e + v e v^u^) 
go into neutrinos. 

Taken together, the neutrino flux connected to the ob- 
served cosmic ray spectrum is given by: 

« 2.3 X UT^e,^ GeV cm -2 s^sr" 1 , (3) 

where the parameter £z accounts for the effects of red- 
shift dependent source evolution. Waxman and Bahcall 
originally presented this argument as an upper bound, 
derived for the case of = 1 (representing the maximum 
flux from a class of optically thin sources). 1 At produc- 
tion, the neutrino flux from positively charged pion de- 
cays consists of equal fractions of v ei v ^ and T>^. After 
oscillations are taken in account, however, the muon neu- 
trinos and anti-neutrinos become a roughly equal mixture 
of muon and tau flavors. 

It should be noted that if we drop the assumption 
that the cosmic ray injection spectrum is of the form, 
dNcR/dEcn oc Bcr.' we cou ld potentially increase the 
flux of 1-10 PeV cosmic neutrinos, although not by more 
than a factor of about two [B] . More significantly, if the 
sources of the ultra high-energy cosmic rays are not opti- 
cally thin to photo-meson production, as assumed in the 
derivation of the Waxman-Bahcall bound, then the ob- 
served cosmic ray spectrum could be made up of only the 
tail of the distribution of accelerated cosmic rays which 
escape from their sources. In such a case, much or even 
most of the energy that goes into accelerating ultra high- 
energy cosmic rays could be lost to the source environ- 
ment, reabsorbing most would-be cosmic rays before they 



1 This result was derived specifically for photo-meson interactions 
and would be increased by a factor of 4/3 if pp interactions had 
instead been assumed. 
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escape. From a class of such optically thick "hidden" 
sources, neutrino fluxes could plausibly be produced that 
are in excess of the Waxman-Bahcall bound by more than 
an order of magnitude, without exceeding the observed 
cosmic ray spectrum or gamma-ray background [BJ. 

So far, this calculation has assumed that the cosmic 
ray spectrum consists of only protons, rather than iron 
or other nuclei species. Heavy nuclei cosmic rays only 
produce charged pions after disintegrating into their con- 
stituent nucleons. Thus if the cosmic ray spectrum were 
dominated by heavy nuclei, the resulting neutrino flux 
would be further suppressed relative to the maximum 
value found in the standard Waxman-Bahcall calcula- 
tion [TJ. In the energy range around ~10 17 eV that leads 
to the production of ~1-10 PeV neutrinos, however, there 
is evidence that the cosmic ray spectrum is dominated by 
protons (although this appears to gradually evolve to a 
heavy nuclei dominated spectrum above 10 19 eV) [HUH]- 
We thus assume a proton dominated cosmic ray spectrum 
throughout this study. 

For a high-energy cosmic neutrino spectrum, we can 
calculate the estimated rate of PeV-scale showers ex- 
pected at IceCube. Hadronic showers can be generated 
through the neutral current interactions of all neutrino 
flavors, with a typical shower energy that is about a quar- 
ter of that possessed by the initial neutrino. Alterna- 
tively, charged-current interactions of electron neutrinos 
and anti-neutrinos produce a superposition of electro- 
magnetic and hadronic showers which collectively con- 
tain the entire energy of the incoming neutrino. While 
electromagnetic and hadronic showers are, in principle, 
distinguishable by their respective muon content, such a 
separation is generally expected to be difficult. In addi- 
tion to producing tracks associated with charged leptons, 
charged current interactions of muon and tau neutrinos 
produce hadronic showers similar to those resulting from 
neutral current processes. 

PeV showers appear to the IceCube detector as photo- 
electrons distributed over an approximately ^300 m ra- 
dius sphere. Although the two recently reported shower 
events were entirely contained within the volume of the 
experiment, IceCube should be capable of detecting par- 
tially contained showers as well. 

The probability that a neutrino passing through the 
effective area of IceCube produces an observable shower 
via a neutral current interaction is given by 
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where a is the neutrino- nucleon cross section [10] . y is 
the energy fraction transferred from the initial neutrino 
to the target nucleus, and L is the length of the detector. 
For charged-current electron neutrino interactions, all of 
the neutrino energy goes into the shower, leading instead 

tO P,y^showor — pNACrL. 

In calculating the expected rate of PeV-scale events at 
IceCube, we must also account for neutrino absorption 
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FIG. 1: To account for the two fully contained PeV shower 
reported by IceCube, a neutrino (plus anti-neutrino, all fla- 
vors) flux roughly at the level depicted by the dotted line is re- 
quired. For comparison, we show (as solid lines) the Waxman- 
Bahcall upper limit on the diffuse neutrino flux |5 , as derived 
for sources which follow a star formation rate-like redshift dis- 
tribution, and for no redshift dependent source evolution. As 
a dashed line, we also show the less stringent limit derived for 
optically thick sources [6]. The error bars represent the spec- 
trum of atmospheric neutrinos, as measured by IceCube [11] . 



in the Earth. A ~1 PeV neutrino has only a few percent 
chance of passing through the equatorial diameter of the 
Earth without undergoing at least one interaction. In 
contrast, for neutrinos with an inclined trajectory of 10 
degrees below the horizon, only about half of such parti- 
cles will undergo one or more interactions in the Earth. 
The effects of absorption are negligible for downward- 
going neutrinos in the energy range under consideration. 
In calculating the effects of neutrino absorption, we adopt 
a simple density model of the Earth, consisting of a 2500 
km radius core with a uniform density of 11,000 kg/m 3 
and a uniform density outer region normalized to the 
overall mass and radius of the Earth. 

For a high-energy cosmic neutrino spectrum that satu- 
rates the Waxman-Bahcall bound (e ff = 1), from sources 
distributed without significant redshift evolution (£z = 
1), we calculate that IceCube should observe a rate of 
13.6 showers per year, per cubic kilometer, with more 
than 1 PeV of energy. If we further require that the 
showers be initiated in the inner volume, restricting our- 
selves to fully contained events, the effective volume of 
IceCube is reduced to roughly ~0.1 km 3 , and thus we 
predict IceCube to observe 1 or 2 fully contained PeV- 
showers per year. 

In Fig. [T] we compare the approximate neutrino flux 
required to explain the two PeV showers reported by 
the IceCube collaboration to the Waxman and Bahcall 
bound, both without source evolution and with source 
evolution that follows the star formation rate. We also 
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show the less stringent upper bound derived in the case 
of optically thick sources [6]. The fact that the neutrino 
flux required to explain IceCube's two PeV showers lies 
not far below the Waxman-Bahcall bound implies that a 
significant fraction (£10%) of energy injected into ~10 17 
eV cosmic rays must be transferred into pions. We also 
note that IceCube's current (90% CL) upper limit on the 
diffuse neutrino flux in this energy range is 2.7 x 10~ 8 
GeV cm~ 2 s _1 sr _1 , assuming an E~ 2 spectrum between 
3.5 x 10 4 and 7 x 10 6 GeV [12], which is very close to the 
flux required to account for the two reported events. 

In the following section, we move away from these gen- 
eral arguments, and consider gamma-ray bursts as a spe- 
cific class of sources that is particularly well suited to 
produce the required flux of neutrinos in the PeV energy 
range. 



III. PEV NEUTRINOS FROM GAMMA-RAY 
BURSTS 

Gamma-ray busts (GRBs) constitute one of the most 
promising sources of high and ultra-high energy cosmic 
rays, and may be capable of accelerating protons to ener- 
gies as high as ~10 20 eV [13] ■ Furthermore, as their name 
implies, gamma-ray burst fireballs contain high densities 
of gamma-rays, enabling for the efficient production of 
neutrinos via the photo-meson interactions of high en- 
ergy protons [HI [TS] . 

More specifically, typical GRBs exhibit a broken 
power-law spectrum of the form: dN^/dE^ cx E~ for 
E 1 £0.1-1 MeV and dN 7 /dE 7 cx E~ x at lower ener- 
gies |16j . Furthermore, the radiation pressure resulting 
from the very high optical depth of GRB fireballs leads to 
their ultra-relativistic expansion, accelerating the plasma 
to Lorentz factors on the order of T ~ 10 2 — 10 3 . In order 
for proton-photon collisions in this environment to ex- 
ceed the threshold for pion production, the proton must 
have an energy (in the observer's frame) that meets the 
following condition (TTl [18] : 
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where z is the redshift of the burst. For any falling spec- 
trum of high-energy protons, such interactions will pre- 
dominantly take place near this threshold. After taking 
into account that only about 20% of the proton's energy 
goes into the charged pion produced in such an interac- 
tion, and that each neutrino carries away only about a 
quarter of the charged pion's energy, this leads to the 
production of neutrinos of characteristic energy: 
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Thus for protons interacting with photons near the ob- 
served spectral break, the resulting neutrinos will have 



energies near that of the two events reported by IceCube. 
The neutrino flux at energies below this value will be 
suppressed by the lack of sufficiently high energy target 
photons in the fireball. For this reason, the PeV energy 
scale is where one roughly expects to observe the first 
GRB neutrinos. 

The overall normalization of the diffuse neutrino flux 
from all GRBs depends on how much of the bursts' in- 
ternal energy goes into accelerating protons to energies 
of ~10 16 eV and above. It has been appreciated for 
some time that if the majority of the highest energy 
cosmic rays originate from GRBs, then one should also 
expect an observable diffuse flux of high energy neutri- 
nos [T7JHH]- In our calculations, we adopt a ratio of ten 
between the energy that goes into accelerated protons 
and electrons (THIHO]; and also assume that 10% of the 
energy in accelerated protons goes into neutrinos, corre- 
sponding to e n ~ 0.2. The actual value of could easily 
be larger or smaller than this value by a factor of several, 
and likely varies significantly from burst-to-burst [15] , 

Most observed GRBs exhibit maximum isotropic lu- 
minosities in the range of L max ~ 10 51 — 10 53 erg/s. 2 
Such bursts are referred to in literature as high luminos- 
ity GRBs and are further divided into short and long du- 
ration bursts with timescales of 0.1-1 and 10-100 seconds, 
respectively; with the majority of observed bursts being 
of long duration [22 24J. In addition, a few low luminos- 
ity GRBs with L max ~ 10 47 erg/s have been detected, 
potentially representing another distinct GRB popula- 
tion [35] (see also Ref. [H]). These low luminosity GRBs, 
which are potentially much more numerous than their 
high luminosity counterparts, generally exhibit smooth 
light curves, wider emission cones and longer durations 

(50-1000 s) pans]. 

Both high and low luminosity GRB populations can be 
described by a luminosity distribution parametrized as: 



(7) 



where <E>o normalizes the luminosity distribution to 1 
(within 2 orders of magnitude above/below the luminos- 
ity break Lb). At present, there is considerable varia- 
tion in the values of these parameters as they appear in 
the literature [25[ [371 US] • In our calculations, we adopt 
ranges based on 2 years of Swift data as presented in 



2 As the observed hard x-ray and 7-ray luminosity is synchrotron 
emission from internal shocks in the relativistic fireball 1211 . this 
emission will be relativistically beamed to within an opening an- 
gle on the order of 9 ~ 1/T. In our calculations, we use the 
isotropic equivalent luminosity related to the true luminosity by: 
Li so = -Ltrue/(1 — cos0). This avoids overestimating the neu- 
trino or photon fluxes. We also estimate the total isotropic en- 



ergy emitted to be E iso ~ Z/™* x T dur , with r du 



2 sec for high 



luminosity GRBs and and r dm . = 50 sec for the low luminosity 
sample. Here, T dur represents duration between which 25% and 
75% of the total energy has been emitted. 



5 



Ref. [25]: L b = (1.2 ± 0.6) x 10 52 erg/s, a x = 0.65 ±0.15, 
a 2 = 2.3 ± 0.3 for high luminosity GRBs and Li, — 
(1.0 ± 0.3) x 10 47 erg/s, a x = ± 0.5, a 2 = 3.5 ± 0.5 for 
low luminosity GRBs. These parameters lead to mean 
luminosities of L = 3.7 x 10 51 to 3.8 x 10 52 erg/s for high 
luminosity GRBs and L = 3.5 x 10 46 to 1.6 x 10 47 erg/s 
for low luminosity GRBs. Throughout the remainder of 
this study, we will use as reference values L^ x = 1 x 10 51 
erg/s and = 10 47 erg/s, but will also consider varia- 
tions in order to test the assumptions of Ref. [27] and [25] . 

For the redshift distribution (co-moving rate density) 
of GRBs, we adopt the star formation rate of Ref. [3D] 
(model "SF2"): 
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High luminosity GRBs occur at a local rate on the order 
of 1 Gpc _3 yr _1 [23 [31] , while the rate of low luminosity 
GRBs is considerably higher, with estimate ranging from 
230 [32] to 5 x 10 3 [33] Gpcr 3 yr -X - We take as reference 
values i?GREs(°) = 1-0 and Rg% b (0) = 350 Gpc- 3 yr- 4 . 
While it has been suggested [29] that the GRB rate red- 
shift distribution does not follow the star formation rate, 
but is instead suppressed at z > 3, this would not affect 
the observed fluxes from high luminosity GRBs by more 
than about a factor of 2 from our reference values for 
high luminosity GRBs. 

The diffuse flux of neutrinos or photons at the location 
of the Earth from the population of all GRBs is given by: 
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and dN ohs /dE obs refers to the observable neu- 
trino/photon flux from an individual GRB located at 
distance D(z): 
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We take as z max = 9.4, which is the most distant GRB 
that has been detected [54] , 

The spectrum of neutrinos (and anti-neutrinos) at in- 
jection can be approximated by a doubly broken power- 
law 1351: 
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The first of these spectral features (at E v = E±) cor- 
responds to the pion production threshold for scatter- 
ing off photons at the observed break in the gamma-ray 
spectrum of GRBs, while the higher energy break (at 
E v = E 2 ) appears as a result of the synchrotron cool- 
ing of muons and pions. The exact locations of these 
breaks is different for typical high and low luminosity 
GRBs due to differences in the strengths of the fireballs' 
magnetic and radiation fields. In calculating the loca- 
tion of the first break, we assume a gamma-ray spectrum 
which breaks at 1 McV or 0.1 MeV for high and low 
luminosity GRBs, respectively. 

In Fig. [2] we show the diffuse flux of neutrinos and anti- 
neutrinos from GRBs for our default parameter choices, 
and for some representative variations of these parame- 
ters. The solid and dashed lines represent the predicted 
flux for our default assumptions, from high and low lu- 
minosity bursts, respectively. The dot-dashed line repre- 
sents the contribution from high luminosity GRBs with 
a redshift distribution which is suppressed above z = 3. 
The dotted line shows the flux from high luminosity GRB 
with alternative choices for the parameters leading to the 
location of the spectral breaks. While these variations 
are by no means exhaustive, they demonstrate that for a 
fairly wide range of assumptions, GRBs are expected to 
generate fluxes of PeV neutrinos that are similar to that 
implied by IceCube's two events (see the dotted line in 
Fig. 1). In particular, for each of the four neutrino fluxes 
shown in Fig. [2] one expects a rate of 4-7 showers with 
energies above 1 PeV per cubic kilometer, per year. For 
an estimate of ~0.1 km 3 for the effective volume for fully 
contained showers, this range of rates is within about 
a factor of two of that required to account for the two 
events reported by IceCube. 

Among the various uncertainties that are involved in 
the calculation of the neutrino flux from GRBs, a few 
stand out as particularly important. First of these is the 
fraction of the energy that goes into accelerating protons 
to energies above 10 16 eV. In the case of the contribu- 
tion from high luminosity GRBs, uncertainties in the lu- 
minosity function are particularly significant. For low 
luminosity GRBs, the local rate of such objects consti- 
tutes a major uncertainty. At present, these uncertainties 
can collectively impact expectations for the diffuse neu- 
trino flux by as much as an order of magnitude and will 
continue to do so until they are better observationally 
constrained. 

By restricting an analysis to events which correlate in 
time and/or direction to known GRBs, it is possible to 
conduct a nearly background free search for neutrinos 
originating from GRBs. Recently, the IceCube collabo- 
ration has applied such a strategy, and used the results 
to derive a stringent upper limit on the flux of high en- 
ergy neutrinos from observed GRBs [36] . Under standard 
astrophysical assumptions, this limit implies that GRBs 
cannot be the only sources of the highest energy (> 10 18 
eV) cosmic rays (see also, Ref. [37]). The two events be- 
ing considered in this paper, however, could still originate 
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FIG. 2: The contribution of gamma-ray bursts (GRBs) to the 
diffuse neutrino (plus anti-neutrino) spectrum. Results are 
shown for high luminosity (solid) and low luminosity (dashed) 
GRBs, calculated using our default parameters, and for high 
luminosity GRB models with a suppressed high redshift dis- 
tribution (dot-dash) and alternative spectral characteristics 
(dots). Each of these models yields a rate of PeV events which 
is comparable to that implied by the two events reported by 
IceCube. See text for details. 




E v (GeV) 

FIG. 3: The contribution of active galactic nuclei (AGN) to 
the diffuse neutrino (plus anti-neutrino) spectrum. Results 
are shown for the models of Protheroe [ID] and Stecker et 
al. gl]. 



from GRB if either, 1) a greater fraction than expected 
of the high energy neutrinos from GRBs originate from 
bursts which are not sufficiently luminous to be observed 
by gamma-ray or x-ray observatories, or 2) a significant 
fraction of the 10 16 — 10 18 eV cosmic ray spectrum (but 
not most of the > 10 18 eV cosmic ray spectrum) orig- 
inates from GRBs. This later possibility is attractive 
from the perspective of the cosmic ray spectrum's chem- 
ical composition. Measurements from the Pierre Auger 
Observatory of the depth of shower maxima and its vari- 
ation suggest that the highest energy cosmic rays are 
largely of heavy chemical composition (closer in mass to 
iron nuclei than protons) , while the composition becomes 
steadily lighter at lower energies, appearing to be domi- 
nated by protons at 10 18 eV [55] • As ultra high-energy 
nuclei accelerated in a GRB are expected to be entirely 
disintegrated into individual nucleons before escaping the 
fireball [7] , the possibility that GRBs provide much of the 
cosmic ray spectrum below ~I0 18 eV, but that another 
class of sources provide the bulk of the highest energy 
(heavy nuclei) cosmic rays, is a well motivated one. 

The fact that the two PeV shower events reported by 
IceCube do not correlate in time with any known GRB 
does not significantly disfavor the hypothesis that these 
events originate from GRBs. The ability of gamma- 
rays and x-ray observatories to monitor for GRBs is, at 
present, substantially incomplete. In particular, the Swift 
Burst Alert Telescope (BAT) and the Fermi Gamma-ray 
Burst Monitor (GBM) collectively cover less than ~2/3 of 
the sky at any given time. With this in mind, one cannot 



rule out the non-negligible possibility that the two PeV 
shower events originated from GRBs which happened to 
fall outside of the combined field-of-view of these instru- 
ments. Furthermore, many GRBs, while in the field-of- 
view of either Swift's BAT or Fermi's GBM, may still go 
undetected if they are of sufficiently low luminosity, or 
are sufficiently distant. 



To assess whether this later possibility is consistent 
with a GRB interpretation of these events, we must at- 
tempt to estimate the efficiency with which neutrino- 
producing GRBs are detected in the gamma-ray and x- 
ray bands 22 , 23 . Given the fluence sensitivity of the the 
Swift's BAT and Fermi's GBM -1 x 1(T 8 erg cm" 2 s" 1 
and ^2 x 10~ 8 erg cm _2 s _1 (at 20 keV) respectively, 
one can estimate how distant a GRB (of a given lumi- 
nosity) could be and still trigger these detectors. We find 
that these experiments should be capable of detecting 
essentially all high luminosity GRBs (L > 10 51 erg/s) 
within their fields-of-view out to a distance of about 8 
Gpc (z « 5). Thus the observed collection of high lumi- 
nosity GRBs is fairly complete (within the given fields- 
of-view). In contrast, low luminosity GRBs (L ~ I0 47 
erg/s) are likely to be detected only within a radius of 
~I00 Mpc, suggesting that the vast majority of the dif- 
fuse neutrino flux from low luminosity GRB will not be 
correlated in time or direction with any observed gamma- 
ray or x-ray signal. 
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IV. PEV NEUTRINOS FROM OTHER 
ASTROPHYSICAL SOURCES 



A. Active Galactic Nuclei 



The kinematics of high-energy neutrino production in 
Active Galactic Nuclei (AGN) is similar to that described 
for GRBs in Eqns. [5][6j but with potentially important 
differences. In particular, the Lorentz factors of AGN 
jets are significantly lower than those of GRB shocks; 
values of T ~ 30 rather than ~300 are thought to be typ- 
ical [35]. As a result, ~10-100 PeV protons can exceed 
the threshold for pion production much more easily, re- 
quiring only the presence of ~keV photons (rather than 
the ~100 keV photons required in GRBs). 

In GRBs, the observed photon spectral break (~0.1-1 
MeV) leads to a break at ~1 PeV in the neutrino spec- 
trum. Thus we may expect the first detections of GRB 
neutrinos to appear at around this energy scale. In con- 
trast, AGN do not typically exhibit a spectral peak at 
keV energies, but instead in the UV, typically at around 
~10 eV. This leads one to expect the neutrino spectrum 
to peak EeV energies, much higher than that from GRBs. 
There is a considerable degree of model dependence in 
this conclusion, however, deriving in large part from un- 
certainties in the spectrum of the target radiation fields. 

We illustrate the nature of this uncertainty in Fig. |3j 
with a comparison of two canonical models of the dif- 
fuse neutrino emission from AGN. First, as a solid line, 
we show the diffuse neutrino spectrum as predicted by 
Protheroe [40 . In this model, the scattering of ultra- 
high energy protons with UV radiation leads to a neu- 
trino spectrum which peaks at EeV energies. For this 
spectral shape, most showers initiated within IceCube's 
volume will be of energy 20 PeV or greater. If IceCube's 
existing data does not contain at least a few enormous 
(non-contained) showers of this energy, this AGN model 
will not be able to account for the two reported PeV 
events. In contrast, the model of Stecker et al. [H] pre- 
dicts a neutrino spectrum from AGN which peaks at a 
much lower energy of a few PeV, not unlike the predic- 
tions from GRBs. This is in large part due to the high 
density of ambient x-rays present in the Stecker et al. 
AGN model. 

We also point out that while neutrino emission from 
known GRBs can be efficiently constrained by searching 
in the time window around the occurrence of a given 
burst, such a background-free strategy is not possible 
for AGN. As a result, while it may be possible to rule 
out (high luminosity) GRB as the source of IceCube's 
two PeV events by searching in the time and direction 
of known GRBs, it will be much more difficult to defini- 
tively test the hypothesis that these neutrinos originate 
from AGN. 



B. Starburst Galaxies 

Collisions of cosmic rays with the radiation in galaxies 
undergoing periods of rapid star formation are predicted 
to yield significant fluxes of ~TeV-PeV neutrinos [42j . 
For reasonable assumptions, and calibrating to the ob- 
served fluxes of gamma-rays [43 and radio emission from 
such objects, contributions to the diffuse neutrino flux are 
expected to be on the order of ~1 x 10~ 8 GeV cm -2 s _1 
sr _1 [42 44 . While it is not clear whether this diffuse 
spectrum extends to energies as high as ~PeV, if it does, 
starburst galaxies could potentially contribute at a level 
capable of accounting for the two PeV events reported 
by IceCube. 

C. Cosmogenic Neutrinos 

In order to generate ~1-10 PeV neutrinos via the pro- 
duction and decay of charged pions, protons with energies 
on the order of of ~10 2 PeV are required. Furthermore, 
in order to exceed the threshold for pion production in 
the center-of-mass frame, these ultra-energetic protons 
must encounter a sufficient density of energetic photons 
(or nucleons). A 10 17 eV proton, for example, can only 
produce pions by scattering off of photons in excess of 
approximately 10 eV. As the background densities of ra- 
diation at such energies are insufficient to lead to any 
significant amount of pion production [2! , we restate the 
well known conclusion that the universe is transparent to 
10 17 eV protons, and thus find that cosmogenic neutrinos 
are unlikely to be the source of the two events reported 
by IceCube. 

We also note that while electron anti-neutrinos from 
the decays of ultra-high energy neutrons could, in prin- 
ciple, contribute to the diffuse neutrino flux at PeV en- 
ergies, such neutrinos are insufficient to account for the 
events reported by IceCube, even in heavy nuclei domi- 
nated cosmic ray scenarios [21 [3] • 

Finally, large scale intergalactic shock waves produced 
due to structure formation may be capable of acceler- 
ating both electrons and protons to Lorentz factors as 
high as ~10 7 [13 US]. Through inverse Compton scat- 
tering, the resulting cosmic ray electrons can produce up 
to 10% of the observed isotropic gamma-ray spectrum 
above 50 GeV. The corresponding protons, however, are 
not predicted to produce significant fluxes of gamma-rays 
or neutrinos. 



V. PREDICTIONS FOR FUTURE 
OBSERVATIONS 

If the two PeV shower events observed by IceCube are, 
in fact, neutrino induced events, then one should expect 
them to be accompanied by a significant rate of corre- 
sponding muon track events. Once produced through 
the charged current interaction of a muon neutrino (or 
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antineutrino), a muon loses energy at a rate given by: 
dE„ 



(13) 



where, in ice, a w 0.002 GeV cm 2 /g and ft « 4.2 x 
10~ 6 cm 2 /g. Thus, before the energy of a propagating 
muon falls below -E' hr , it travels a distance of: 



log 



(14) 



For PeV-scale muons, this range will significantly exceed 
the linear scale of IceCube, increasing the effective tar- 
get mass of the experiment. The probability of a muon 
neutrino producing a muon track with an energy greater 
than i?* hr inside of the volume of IceCube is given by: 



Pu 



> track 



-(E?*/E„ 



y))^j-dy. 
ay 



To make a comparison between the rate of shower and 
muon track events, we consider a spectrum of the form 
dN v /dE v cx E~ 2 , and with the flavor ratios predicted 
from the decays of positively charged mesons, after ac- 
counting for oscillations (y e : fn : v T : D e : : v T ss 
2:1:1:0:1:1). 3 Considering only neutrinos with 
energies in excess of 1 PeV, we find that for every fully 
contained shower (V B g « 0.1 km 3 ), one expects 33 (11) 
muons tracks with more than 100 TeV (1 PeV) of en- 
ergy to travel through the volume of IceCube. We thus 
conclude that if these two fully contained PeV showers 
are neutrino induced, then the existing data set from Ice- 
Cube should contain tens of PeV muon track events, as 
well as tens of partially contained PeV showers. Such 
events would be very valuable for searching for correla- 
tions in direction (for track events) and time with known 
gamma-ray bursts. 

We also note that it may be possible for IceCube to 
identify PeV-scale events which are unique to tau neu- 
trinos [IS]. Following the approach of Ref. [JS], and for 
the same assumptions as in the previous paragraph, we 
estimate that for every fully contained PeV shower, one 
should expect approximately 0.7 lollipop events and 0.25 
double bang events. 4 While these are not particularly 
high rates, even one or two tau-unique events would be 
very helpful in confirming an astrophysical nature of the 
events in question. 

VI. SUMMARY AND CONCLUSIONS 

Experiments such as IceCube were motivated in large 
part by the connection between the observed cosmic ray 



spectrum and high-energy cosmic neutrinos. In partic- 
ular, if a significant fraction of high-energy cosmic rays 
interact via photo-meson interactions, then it should be 
possible to observe the resulting diffuse flux of very high- 
energy cosmic neutrinos. With the observation of two 
showers with energies in excess of 1 PeV at IceCube, we 
may very well be witnessing the first light in the field of 
high-energy neutrino astronomy. And if so, the energy 
and flux implied by these events is strongly suggestive of 
a connection with the observed cosmic ray spectrum. 

In this paper, we have considered a number of pos- 
sible sources and mechanisms that could produce PeV 
neutrinos in a quantity capable of accounting for the two 
events observed by IceCube. We find that photo-meson 
interactions of ~10-100 PeV protons to be the most 
promising possibility, which could plausibly be realized in 
a variety of astrophysical sources, including gamma-ray 
bursts (GRBs), active galactic nuclei (AGN), and star- 
burst galaxies. Due to the observed break at 0.1-1 MeV 
in the spectra of typical GRBs, one expects the result- 
ing neutrino spectrum from such objects to peak at the 
PeV scale; in fact we find that GRBs can naturally ex- 
plain the two observed events. More specifically, we find 
that a neutrino flux comparable to that implied by the 
two events is predicted for a wide range of assumptions 
regarding the redshift distribution, luminosity function, 
and other physical characteristics of GRBs. In contrast, 
the photon spectra in typical AGN models peak at much 
lower energies, and thus neutrino spectra are predicted 
which peak at energies well above those of IceCube's two 
events. In fact, searches for even higher energy showers in 
IceCube's existing data set may be capable of excluding 
some of such models. Cosmic ray collisions in starburst 
galaxies are also a promising possibility, but it is not clear 
whether neutrinos with energies as high as 1-2 PeV can 
be generated in such sources. 

Lastly, we point out that the two fully contained 
shower events observed by IceCube should be accompa- 
nied by tens of PeV muon tracks and partially contained 
showers. A more modest number of events unique to tau 
neutrinos should also be expected to appear. Thus as 
IceCube's overall data set is further scrutinized, it should 
become increasingly possible to determine whether these 
two events are in fact cosmic in nature and, if so, to begin 
to constrain their origin. 
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3 Electron anti-neutrinos will be produced through oscillations of 
Vp, and v T for non-zero values of $13. Thus the rate of showers 
at the energy of the Glashow resonance (6.3 PeV) could be used 
as a probe of this oscillation parameter |47) . 

4 Here, we have required only that the showers involved in these 



events are initiated within the volume of IceCube, and are not 
necessarily fully contained. 
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